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ABSTRACT 

We study the intergalactic transmission of radiation in the vicinity of the Lya wavelength. Simu- 
lating sightlines through the intergalactic medium (IGM) in detailed cosmological hydrosimulations, 
the impact of the IGM on the shape of the line profile from Lya emitting galaxies at redshifts 2.5 to 
6.5 is investigated. In particular we show that taking into account the correlation of the density and 
velocity fields of the IGM with the galaxies, the blue part of the spectrum may be appreciably reduced, 
even at relatively low redshifts. This may in some cases provide an alternative to the often-invoked 
outflow scenario, although it is concluded that this model is still a plausible explanation of the many 
asymmetric Lya profiles observed. 

Applying the calculated wavelength dependent transmission to simulated spectra from Lya emitting 
galaxies, we derive the fraction of photons that are lost in the IGM, in addition to what is absorbed 
internally in the galaxies due to dust. 

Moreover, by comparing the calculated transmission of radiation blueward of the Lya line with 
corresponding observations, we are able to constrain the epoch when the Universe was reionized to 
z < 8.5. 

Subject headings: intergalactic medium — radiative transfer — scattering — line: profiles 



1. INTRODUCTION 

The past decade has seen a rapid increase in the use 
of Lya as a cosmological tool. A plethora of physical 
characteristics of young galaxies and the intergalactic 
medium (IGM) that separate them manifest themselves 
in the shape, the strength, and the spatial distribution 
of both Lya emission and absorption lines. 

Since generally the evolution of astrophysical objects 
happen on timescales much longer than a human life- 
time, numerical simulations are particularly well-suited 
for exploring the Universe. Many elaborate numerical 
codes have been constructed to model the Universe on 
all scales, from dust agglomeration, over planetary and 
stellar formation, to simulations of cosmological volumes. 
Although these codes include a steadily increasing num- 
ber of physical processes, numerical resolution, etc., and 
although they are able to predict various observables, 
many of them fail to account for the fact that the ob- 
served light may be rather different from the light that 
was emitted. Once radiation leaves its origins, it may 
still be subject to various physical processes altering not 
only its intensity, but in some cases also its spatial and 
spectral distribution. If we do not understand these pro- 
cesses, we may severely misinterpret the predictions of 
the models when comparing to observations. This is 
particularly true with regard to resonant lines like Lya, 
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where its complicated path may take it through regions 
of unknown physical conditions. In order to interpret the 
observed Lya lines correctly, it is crucial to have an un- 
derstanding of the physical processes that influence the 
shape and transmission of the line. 

Although the atomic processes regulating the individ- 
ual scatterings are always the same 6 , the physical condi- 
tions governing different regions in space will have a large 
impact on the outcome. For example, photons produced 
by gravitational cooling are mostly born in low-density 
regions, and may escape the galaxy more or less freely. 
These photons thus tend to have a frequency close to the 
line center. On the other hand, photons produced in the 
dense stellar regions will usually have to scatter far from 
the line center in order to escape, and thus tend to com- 
prise the wings of the spectrum. However, these photons 
are also more vulnerable to dust, since their paths are 
much longer, and since also most of the galactic dust 
resides in these emission regions. Moreover, the veloc- 
ity field of the gas elements also has an effect on the 
observed line profile. In particular macroscopic gas mo- 
tions, such as gas accretion or outflows caused by feed- 
back due to massive star formation, will skew the line 
(e.g. Dijkstra et al. 2006; Verhamme et al. 2006). 

While the history of theoretical Lya radiative trans- 
fer (RT) dates back to the 1930's (Ambarzumian 1932; 
Chandrasekhar 1935), the first practical prediction of 
spectra emerging from a source embedded in neutral hy- 
drogen was given by Harrington (1973), and later gen- 
eralized by Neufeld (1990). Due to the high opacity 
for a photon at the line center, in general the photons 
have to diffuse in frequency and should hence escape the 
medium in a double-peaked spectrum. Although this 

6 As long the atom is not perturbed during the ~10~ 8 seconds it 
is excited; this, however, can be safely ignored in most astrophysical 
situations. 
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has indeed been observed (e.g. Yee & De Robertis 1991; 
Venemans et al. 2005; Vanzella 2008), apparently most 
Lya profiles from high-rcdshift galaxies seem to be miss- 
ing the blue peak. An immediate conclusion would be 
that high-redshift, Lya emitting galaxies (LAEs) are in 
the process of massive star formation and thus exhibit 
strong outflows. Indeed, this scenario has been invoked 
to explain a large number of LAE spectra, most convinc- 
ingly by Verhamme et al. (2008) who, assuming a central 
source and thin surrounding shell of neutral gas while 
varying its expansion velocity, temperature, and gas and 
dust column density, manage to produce nice fits to a 
number of observed spectra. 

Although the thin expanding shell scenario hinges on 
a physically plausible mechanism, it is obviously rather 
idealized. Furthermore, since most observations show 
only the red peak of the profile, it seems to indicate that 
most (high-redshift) galaxies exhibit outflows. However, 
at high redshifts many galaxies are still forming, resulting 
in infall which would in turn imply an increased blue 
peak. Since this is rarely observed, in this paper we aim 
to investigate whether some other mechanism, such as 
IGM RT effects, could be responsible for removing the 
blue peak and/or enhancing the red peak. 

The present study complements other recent endeav- 
ors to achieve a comprehensive understanding of how 
the Lya line is redistributed in frequency and real space 
after having escaped its host galaxy (Iliev et al. 2008; 
Barnes & Haehnelt 2010; Faucher-Giguere et al. 2010; 
Zheng et al. 2010a,b,c). Our cosmological volume is not 
as large as most of these studies, and although gas dy- 
namics are included in the simulations, ioniziation by UV 
radiation is calculated as a post-process rather than on 
the fly. The reward is a highly increased resolution, al- 
lowing us to study the circumgalactic environs in great 
detail. Moreover, we inquire into the temporal evolution 
of the IGM. 

The exploration of the high-redshift Lya Universe is 
rapidly progressing, gradually providing a census of the 
physical properties of young galaxies, such as their lu- 
minosity functions, clustering properties, evolution, and 
contents. Recent notable surveys have revealed hundreds 
of LAEs at high redshifts, from z ~ 2.1-2.3 (Guaita et al. 
2010; Nilsson et al. 2009), to z ~ 4.5 (Wang et al. 2009), 
to z ~ 6.6 (Ouchi et al. 2010), and for small samples 
even farther (Hibon et al. 2010; Tilvi et al. 2010), and 
the present work, with its emphasis on the reshaping of 
the Lya line profile and amplitude throughout different 
epochs is a step toward a more accurate interpretation 
of the observations. 

In the following section, we discuss qualitatively the 
impact of the IGM on the Lya line profile. In Section 3 
we then decribe the numerical simulations we have per- 
formed to study the effect quantitatively, the results of 
which are presented in Section 4. We discuss our findings 
in Section 5, and summarize in Section 6. 

2. THE INTERGALACTIC MEDIUM 

As light travels through the expanding Universe, it gets 
rcdshiftcd, implying that wavelengths A blucward of the 
Lya line center are eventually shifted into resonance. If, 
for a given wavelength, this happens in the vicinity of a 
sufficient amount of neutral hydrogen, the spectrum ex- 
periences an absorption line (although strictly speaking 



the photons are not absorbed, but rather scattered out 
of the line of sight). This results in the so-called Lya 
forest (LAF; Lynds 1971; Sargent et al. 1980). 

As the red part of the spectrum is only shifted far- 
ther away from resonance, IGM absorption tends to 
skew the line and not simply diminish it by some fac- 
tor. In most earlier studies of the line profiles of high- 
rcdshift galaxies, the IGM has either been ignored (e.g. 
Verhamme et al. 2006, 2008), taken to transmit the red 
half and remove the blue part (e.g. Finkelstein et al. 
2008), or to influence the line uniformly by a fac- 
tor e~^ T \ where (r) is the average optical depth of 
the Universe (e.g. Bruscoli et al. 2003; Meiksin 2005; 
Faucher-Giguere et al. 2008; Rykoff et al. 2009). 

Consider a source emitting the normalized spectrum 
-Fem(A) = 1 for all A. In an idealized, completely ho- 
mogeneous universe undergoing completely homologous 
expansion ( "ideal Hubble flow" ) and with the absorption 
profiles having a negligible width, the observed spectrum 
would simply be a step function, with F(X) = F TC d = 1 
for A > Ao, and F(X) = Fbi uc < 1 for A < Ao, where 
Ao = 1216 A is the central wavelength of the line. Three 
factors contribute to make F differ from a step function: 

Firstly, the Lya line is not a delta function, but has a fi- 
nite width, in some cases extending significantly into the 
red part of the spectrum. Secondly, the IGM is highly in- 
homogeneous, leading to large variations in F blueward 
of Ao — the LAF. If we consider the average effect of 
the IGM and let F be the average of many sightlincs, 
-Fbiuc should still be a constant function of wavelength. 
However, in the proximity of galaxies the gas density is 
higher than far from the galaxies; on the other hand, in 
these regions the stellar ionizing UV radiation may re- 
duce the neutral gas fraction. Consequently, wavelengths 
just blueward of the line center may not on average be 
subject to the same absorption as farther away from the 
line. Regardless of which of the two effects gas over- 
density or stronger UV field — is more important, the 
correlation of the IGM with the source cannot be ne- 
glected. 

Finally the expansion is not exactly homologous, since 
peculiar velocities of the gas elements will cause fluctua- 
tions around the pure Hubble flow. Considering again 
not individual sightlincs but the average effect of the 
IGM, these fluctuations are random and cancel out; on 
average, for every gas element that recedes from a galaxy 
faster than the Hubble flow and thus causes an absorp- 
tion line at a slightly bluer wavelength, another gas el- 
ement does the opposite. Hence, the average transmis- 
sion in a "realistic" universe is the same as in a universe 
where there are no peculiar velocities. However, in the 
proximity of overdensities, the extra mass results in a 
"retarded" expansion of the local IGM. When expansion 
around a source is somewhat slower than that of the rest 
of the Universe, on average matter in a larger region will 
be capable of causing absorption, since the slower the 
expansion, the farther the photons will have to travel 
before shifting out of resonance. 

2.1. The transmission function 

This paper focuses on two different aspects of the in- 
tergalactic transmission. First, to see how the Lya line is 
affected by the IGM, we calculate a "transmission func- 
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tion" given by the average, normalized flux F(X). Specif- 
ically, F is calculated by taking the median value in each 
wavelength bin of many sightlines, originating just out- 
side a large number of galaxies (where "just outside" will 
be defined later). The standard deviation is defined by 
the 16 and 84 percentiles. From the above discussion 
we may expect that F be characterized by a red part 
F rec j — lj an d a blue part F\>\ uc < 1, but with a non- 
trivial shape just blueward of the line center, since the 
IGM in the vicinity of the sources is different than far 
from the sources. 

2.2. The average transmission 

We will also examine the average transmission T = 
T(z) of the IGM, defined by first calculating for each in- 
dividual sightline the fraction of photons that are trans- 
mitted through the IGM in a relatively large wavelength 
interval well away from the line center, and then taking 
the median of all sightlines (again with the 16 and 84 
percentiles defining the standard deviation). This quan- 
tity is sensitive to the overall ionization state of the IGM, 
and has therefore been used observationally to put con- 
straints on the so-called Epoch of Reionization (EoR; e.g. 
Becker et al. 2001; Djorgovski et al. 2001). A paramount 
puzzle in modern cosmology is the question of when and 
how the hydrogen, and later helium, was reionized. The 
EoR marks a comprehensive change of the physical state 
of the gaseous Universe, and to understand the cause, as 
well as the course, of this phenomenon is a challenging 
task. Besides being a compelling event in itself, it also 
has profound implications for the interpretation of ob- 
servations and theoretical cosmological models, not only 
due to the increased transparency of the IGM, but also 
because of the accompanying rise in IGM temperature. 

One way of probing the EoR is by looking at the spec- 
tra of high-redshift quasars. The intrinsic spectrum of a 
quasar is characterized by broad emission lines; in partic- 
ular the Lya line is very prominent. For larger redshifts 
the absorption lines of the LAF gets increasingly copious, 
until eventually all of the light blueward of Lya is ab- 
sorbed, resulting in the so-called Gunn-Petcrson trough 
(Gunn & Peterson 1965). Observations of a large num- 
ber of quasars show that the Universe was largely opaque 
to radiation blueward of Lya at z > 6 (Songaila 2004; 
Fan et al. 2006). 

3. SIMULATIONS 

Laursen et al. (2009a,b) considered numerically the 
shape of the Lya spectrum emerging from galaxies at 
z ~ 3.5, based on hydro/gravity galaxy simulations of 
very high resolution. In these studies, the RT was carried 
out with Monte Carlo simulations, following individual 
photons (or photon packets) as they scatter stochasti- 
cally out through the interstellar medium. Since high 
resolution is important for the Lya RT, especially when 
taking into account dust, the simulations were conducted 
in galaxies extracted from a fully cosmological, large vol- 
ume model, resimulated 7 at high resolution, and inter- 
polated onto an adaptively refined grid. This Lya RT 

7 A "resimulation" is performed by tracing the position of the 
particles comprising a given galaxy back to the initial conditions, 
and then repeating the simulation at 8 or 64 times the original 
mass resolution. 



evidently only predicts the spectrum of radiation that 
one would observe if located in the vicinity of the galax- 
ies. In reality the radiation has to travel trough the IGM. 
At a rcdshift of ~3.5, the IGM is largely ionized, and the 
spectra escaping the galaxies are, in general, not very 
different from what would be observed at Earth. How- 
ever, even a very small amount of neutral hydrogen may 
influence the observations. To investigate just how large 
an impact the IGM exerts on the radiation, full IGM 
RT has to be computed. Furthermore, at higher red- 
shifts where the IGM is generally more neutral and more 
dense, neglecting the effect would lead to severely erro- 
neous results. 

In principle this could be achieved by performing first 
the "galactic" RT in the high-resolution resimulations 
and subsequently continuing the RT in the low-resolution 
cosmological volume from the location of the individual 
galaxies. However, although the physics of scattering 
in galaxies and that of scattering in the IGM is not in- 
herently different, the difference in physical conditions 
imposes a natural division of the two schemes: in the 
dense gas of galaxies, photons are continuously scattered 
in and out of the line of sight, whereas in the IGM, once 
a photon is scattered out of the line of sight, it is "lost" , 
becoming part of the background radiation. The prob- 
ability of a background photon being scattered into the 
line of sight, on the other hand, is vanishingly small. 

In order to disentangle galactic from intergalactic ef- 
fects, and, more importantly, to investigate the general 
effect of the IGM instead of merely the IGM lying be- 
tween us and a handful of resimulated galaxies, we take 
a different approach: the transmission properties of the 
IGM are studied by calculating the normalized spectrum 
F(X) — the transmission function — in the vicinity of 
the Lya line, as an average of a large number of sightlines 
cast through the cosmological volume, and originating 
just outside a large number of galaxies. 

3.1. Cosmological simulations 

A standard, flat ACDM is assumed, with Q m = 0.3. 
n A = 0.7, and h = Ho/100 km - 1 Mpc" 1 = 0.7. 

The radiative transfer is carried out in a cosmo- 
logical volume modeled using self-consistent, ab initio 
hydro/gravity simulations (Sommer-Larsen et al. 2003; 
Sommer-Larsen 2006). The code used in these simula- 
tions is a significantly improved version of the TreeSPH 
code, which has been used previously for galaxy for- 
mation simulations (Sommer-Larsen et al. 2003). The 
main improvements over the previous version are: (1) 
The "conservative" entropy equation solving scheme 
(Springel & Hernquist 2002) has been implemented, im- 
proving shock resolution in SPH simulations. In partic- 
ular, this improves the resolution of the starburst-driven 
galactic supcrwinds invoked in the simulations. (2) Non- 
instantaneous gas recycling and chemical evolution, trac- 
ing the ten elements H, He, C, N, O, Mg, Si, S, Ca and 
Fe (Lia et al. 2002a,b); the algorithm invokes effects of 
supcrnovac of type II and type la, and mass loss from 
stars of all masses. (3) Atomic radiative cooling de- 
pending both on the metal abundance of the gas and 
on the meta-galactic UV background (UVB), modeled 
after Haardt & Madau (1996), is invoked, as well as a 
simplified treatment of radiative transfer, switching off 
the UVB where the gas becomes optically thick to Ly- 
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man limit photons on a scale of 0.1 kpc (Sommer-Larsen 
2006). The simple scheme for UVB RT is supplemented 
by a more elaborate post-process ionizing UV RT scheme, 
as described in Section 3.2. 

The simulation volumes were selected from two dark 
matter (DM) only cosmological simulations of box length 
10/i _1 Mpc, set up at z 2 ; = 39 using identical Fourier 
modes and phases, but with rms linear fluctuations on 
a scale of 8ft. _1 Mpc, erg, equal to either 0.74 or 0.9, re- 
spectively, bracketing the latest WMAP 8 -inferred value 
of -0.8 (Jarosik et al. 2010). 

The original DM-only simulations were run with peri- 
odic boundary conditions using 128 3 DM particles, and 
have been used previously for galaxy resimulations (e.g. 
Sommer-Larsen et al. 2003; Sommer-Larsen 2006, note, 
however, that in these simulations, erg — 1.0 was as- 
sumed). For each of the two cosmological DM-only simu- 
lations, mass and force resolution was increased (by fac- 
tors of eight and two, respectively) in a spherical La- 
grangian sub- volume of radius 5/i _1 Mpc (about half the 
volume of the original DM-only simulation), and in this 
region all DM particles were split into a DM particle and 
a gas (SPH) particle according to an adopted universal 
baryon fraction of fa = 0.15, in line with recent esti- 
mates. The cosmological regions resimulated thus had a 
comoving volume of about 1500 Mpc 3 , contained about 
17 million particles in total, and were resimulated with 
an open boundary condition using the full hydro/gravity 
code. The numerical resolution of these cosmological res- 
imulations was hence the same as that of the individual 
galaxy formation simulations of Sommer-Larsen et al. 
(2003), reflecting the very considerable increase in com- 
puting power since then. 

The masses of SPH, star and DM particles were m gas = 
= 7.3 x 10 s and m DM = 4.2 x 10 6 h~ x M , and 
the gravity softening lengths were e gas = e* = 380 and 
£dm = 680 h^ 1 pc. The gravity softening lengths were 
kept constant in physical coordinates from z = 6 to z = 
0, and constant in comoving coordinates at earlier times. 

The choice of stellar initial mass function (IMF) af- 
fects the chemical evolution and the energetics of stel- 
lar feedback processes, in particular starburst-driven 
galactic supcrwinds. Motivated by the findings of 
Sommer-Larsen & Fynbo (2008), a Salpeter (1955) IMF 
was adopted for the simulations presented in this paper. 

3.2. Ionizing UV radiative transfer 

The Haardt & Madau UVB field becomes negligible at 
z > z le = 6. To comply with the results of WMAP, which 
indicate a somewhat earlier onset of the ionizing back- 
ground, we ran an additional set of models, in which 
the UVB intensity curve was "stretched" to initiate at 
z rc = 10. In the following, Model 1, 2, and 3 refers 
to simulations using (z ro ,er 8 ) = (10,0.74), (10,0.9), and 
(6,0.74), respectively. The models (also summarized in 
Table 1) with z IC = 10 will be referred to as "early" 
reionization, while z TC = 6 will be referred to as "late" 
rcionization. In Sections 4.1 and 4.3 we study the differ- 
ences between these models. 

A realistic UV RT scheme 

(Razoumov & Sommer-Larsen 2006, 2007) was em- 

8 http : //lambda . gsf c . nasa . gov/pr oduct/map/curr ent / 



ployed to post-process the results of the hydrodynamical 
simulations using the following approach: First, the 
physical properties of the SPH particles (density, 
temperature, etc.) are interpolated onto an adaptively 
refined grid of base resolution 128 3 , with dense cells 
recursively subdivided into eight cells, until no cell 
contains more than ten particles. Around each stellar 
source, a number of radial rays arc constructed and 
followed through the nested grid, accumulating photore- 
action number and energy rates in each cell. Rays are 
split either as one moves farther away from the source, 
or as a refined cell is entered. These reaction rates are 
used to compute iteratively the equilibrium ionization 
state of each cell. 

In addition to stellar photons, we also account 
for ionization and heating by LyC photons origi- 
nating outside the computational volume with the 
FTTE scheme (Razoumov & Cardall 2005) assuming the 
Haardt & Madau UVB, modified to match the particular 
rcionization model. 

Since the relative ratio of Hi, Hei, and Hen densities 
varies from cell to cell, the UV RT is computed separately 
in three frequency bands, [13.6, 24.6[ eV, [24.6, 54.4[ eV, 
and [54.4, oo[ eV, respectively. In each cell the angle- 
averaged intensity is added to the chemistry solver to 
compute the ionization equilibrium. 

3.3. Galaxy selection criteria 

Galaxies are located in the simulations as described 
in Sommer-Larsen et al. (2005). To make sure that a 
given identified structure is a real galaxy, the following 
selection criteria are imposed on the sample: 

1. To ensure that a given structure "z" is not just a 
substructure of a larger structure "j" , if the center 
of structure i is situated within the virial radius of 
j, it must have more stars than j. 

2. The minimum number of star particles must be at 
least A^ jm i n = 15. This corresponds to a minimum 
stellar mass of log(M^ jm i n /MQ) = 7.2. 

3. The circular velocity, given by V c = \f GM v - n -/r V i r 
must be V c > 35 km s _1 . 

The IGM in the vicinity of large galaxies is to some ex- 
tend different from the IGM around small galaxies. The 
more mass gives rise to a deeper gravitational poten- 
tial, enhancing the retarded Hubble flow. On the other 
hand, the larger star formation may cause a larger bub- 
ble of ionized gas around it. To investigate the difference 
in transmission, the sample of accepted galaxies is di- 
vided into three subsamples, denoted "small" , "interme- 
diate" , and "large" . To use the same separating criteria 
at all redshifts, instead of separating by mass — which 
increases with time due to merging and accretion — the 
galaxies are separated according to their circular veloc- 
ity, which does not change significantly over time. The 
thresholds arc defined somewhat arbitrarily as V\ = 55 
km s _1 (between small and intermediate galaxies) and 
V2 = 80 km s _1 (between intermediate and large galax- 
ies). The final sample of galaxies is shown in Figure 1, 
which shows the relation between stellar and virial mass, 
and the distribution of circular velocities. 
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Fig. 1. — Left: Scatter plot of virial masses M v - lr vs. stellar masses M* for the full, unfiltcred sample of galaxies. The colors signify 
rcdshift, with more red meaning higher redshift (exact values are seen in the right plot's legend). The data points of Model 1, 2, and 3 
are shown with triangles, plus signs, and crosses, respectively. The galaxies that are rejected are overplottet with gray crosses, with dark, 
medium, and light gray corresponding to rejection criterion 1, 2, and 3, respectively. Right: Distribution of circular velocities V c for the 
accepted galaxies in Model 1. The distributions for Model 2 and 3 look similar, although Model 2 has more galaxies (see Table 1). 



TABLE 1 

Number of galaxies in the simulations 



Model z TC erg z Total Small Intermediate Large 



1. 10 0.74 



2. 10 0.9 



3. 6 0.74 



2.5 


343 


207 


80 


56 


3.5 


309 


138 


127 


11 


3.8 


283 


125 


119 


39 


4.3 


252 


102 


115 


35 


1.8 


225 


111 


86 


28 


5.3 


201 


98 


85 


18 


5.8 


151 


7.5 


62 


17 


6.5 


121 


70 


11 


7 


3.5 


405 


207 


128 


70 


1.3 


384 


165 


150 


69 


1.8 


341 


150 


129 


62 


5.3 


324 


115 


125 


51 


5.8 


277 


121 


106 


50 


6.5 


252 


126 


92 


31 


3.5 


325 


162 


122 


11 


3.8 


318 


165 


117 


36 


1.3 


293 


150 


110 


33 


4.8 


250 


138 


81 


28 


5.3 


204 


101 


85 


18 


5.8 


160 


75 


67 


18 


6.5 


126 


69 


50 
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Note. - "Small", "intermediate", and "large" galaxies are 
defined as having circular velocities V c < 55 km s , 55 km 
s _1 < V c < 80 km s , and V c > 80 km s , respectively. 

The exact number of galaxies in the three models stud- 
ied in this work is seen in Table 1. 

3.4. IGM radiative transfer 

The IGM RT is conducted using the code IGMtrans- 
fer 9 . For the RT we use the same nested grid used for 

9 The code can be downloaded from the following URL: 
http : / /www. dark-cosmology . dk/~pela/IGMtransf er . html . 



the ionizing UV RT. The transmission properties of the 
IGM are studied by calculating the normalized spectrum 
F(X) in the vicinity of the Lya line, as an average of 
a high number of sightlines cast through the simulated 
cosmological volume. 

The resulting value of F(X) at wavelength A for a given 
sightline is 

F(A) = e" T(A) . (1) 

The optical depth r is the sum of contributions from all 
the cells encountered along the line of sight: 

cells 

t(A) = 2J iHi,i n cr(A + \v\\ ti /c). (2) 

i 

Here, n.Hi,i is the density of neutral hydrogen in the i'th 
cell, Ti is the distance covered in that particular cell, v\\j 
is the velocity component of the cell along the line of 
sight, and <r(A) is the cross section of neutral hydrogen. 
Due to the resonant nature of the transition, the largest 
contribution at a given wavelength will arise from the 
cells the velocity of which corresponds to shifting the 
wavelength close to resonance. 

Although no formal definition of the transition from a 
galaxy to the IGM exists, we have to settle on a defini- 
tion of where to begin the sightlines, i.e. the distance r$ 
from the center of a galaxy. Observed Lya profiles result 
from scattering processes in first the galaxy and subse- 
qently the IGM, and regardless of the chosen value of r , 
for consistency "galactic Lya RT" should be terminated 
at the same value when coupling the two RT schemes. In 
galactic Lya RT, individual photons are traced, while in 
IGM RT we only consider the photons that are not scat- 
tered out of the line of sight. Hence the sightlines should 
begin where photons are mainly scattered out of the line 
of sight, and only a small fraction is scattered into the 
line of sight. Since more neutral gas is associated with 
larger galaxies, as well as with higher redshift, clearly 
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ro depends on the galaxy and the epoch, but measuring 
ro in units of the virial radius r v ; r helps to compare the 
physical conditions around different galaxies. In Section 
4.1 we argue the a reasonable value of ro is 1.5r v i r , and 
show that the final results are only mildly sensitive to 
the actual chosen value of ro. 

At a given redshift, F(X) is calculated as the median in 
each wavelength bin of 10 3 sightlines from each galaxy in 
the sample (using 10 4 sightlines produces virtually iden- 
tical results) . The number of galaxies amounts to several 
hundreds, and increases with time (cf. Table 1). The sim- 
ulated volume is a spherical region of comoving diameter 
-Dbox: but to avoid any spuriosities at the edge of the 
region only a 0.9-Dbox sphere is used. Due to the limited 
size of the volume, in order to perform the RT until a 
sufficiently short wavelength is rcdshiftcd into resonance 
the sightlines are allowed to "bounce" within the sphere, 
such that a ray reaching the edge of the sphere is rc- 
"emitted" back in a random angle in such a way that the 
total volume is equally well sampled. Note however that 
in general the wavelength region close to the line that is 
affected by the correlation of the IGM with the source is 
reached well before the the first bouncing. 

The normalized spectrum is emitted at rest wavelength 
in the reference frame of the center of mass of a galaxy, 
which in turn may have a peculiar velocity relative to 
the cell at which it is centered. This spectrum is then 
Lorentz transformed between the reference frames of the 
cells encountered along the line of sight. Since the ex- 
pansion of space is approximately homologous, each cell 
can be perceived as lying in the center of the simulation, 
and hence this bouncing scheme does not introduce any 
bias, apart from reusing the same volume several times 
for a given sightline. However, since the sightlines scatter 
around stochastically and thus pierce a given region from 
various directions, no periodicities arise in the calculated 
spectrum. 

To probe the average transmission, the sightlines are 
propagated until the wavelength 1080 A has been red- 
shifted into resonance, corresponding to Az ~ 0.1. In 
this case the sightlines bounce roughly 30 (20) times at 
z = 2.5 (6.5). 

4. RESULTS 

In the following sections, Model 1 is taken to be the 
"benchmark" model, while the others are discussed in 
Section 5.5. 

4.1. The Lya transmission function 

Figure 2 shows the calculated transmission functions 
F(X) at various rcdshifts. Indeed, a dip just blueward 
of the line center is visible at all redshifts. The results 
in Figure 2 were calculated as the median of sightincs 
emerging from all galaxies in the sample. Of course a 
large scatter exists, since each sightline goes through very 
different regions, even if emanating from the same galaxy. 
Figure 3 shows the scatter associated with F for three 
different redshifts. The equivalent transmission functions 
for the z ro = 6 model are shown in Figure 4. While at 
high rcdshifts the z IC = 6 model clearly results in a much 
more opaque universe, at lower redshifts the transmission 
properties of the IGM in the different models are more 
similar, although at z = 3.5, the z re = 10 model still 
transmits more light. 



To see the difference in transmission around galaxies 
of different sizes, Figure 5 shows transmission function 
for three size ranges (defined in Section 3.3). Since the 
distance from the galaxies at which the sightlines start 
is given in terms of virial radii, sightlines emerging from 
small galaxies start closer to their source than for large 
galaxies, and since at lower redshifts they tend to be clus- 
tered together in the same overdensities as large galaxies, 
this results in slightly more absorption. However, the dif- 
ference is not very significant. 

4.2. Effect on the spectrum and escape fraction 

In Figure 6 the "purpose" of the transmission function 
is illustrated: the left panel shows the spectrum emerging 
from a galaxy of M v ; r = 4.9 x 10 9 M Q and Lya luminosity 
Lhya = 4.9 x 10 40 erg s _1 , calculated with the Monte 
Carlo Lya RT code MoCaLaTA (Laursen et al. 2009a). 
Its circular velocity of 42 km s _1 characterizes it as a 
small galaxy. The spectrum which is actually observed, 
after the light has been transferred through the IGM, 
is shown in the right panel. Although on average the 
effect of the IGM is not very large at this redshift, as 
seen by the solid line in the right panel, due to the large 
dispersion in transmission (visualized by the gray-shaded 
area) at least some such galaxies will be observed with a 
substantially diminished blue peak. 

In general, the larger a galaxy is the broader its emit- 
ted spectrum is, since Lya photons have to diffuse far- 
ther from the line center for higher column densities of 
neutral gas. If dust is present, this will tend to narrow 
the line (Laursen et al. 2009b). Larger galaxies tend to 
have higher metallicities and hence more dust, but the 
lines will still be broader than the ones of small galaxies. 
The galaxy used in Figure 6 is quite small. In Figure 7 
we show the impact of the IGM on the nine simulated 
spectra from Laursen et al. (2009b), spanning a range 
in stellar masses from ~10 7 to ~10 10 . For comparison, 
typical LAEs have stellar masses of M*/Mq ~ 10 8 -10 9 
(Gawiser et al. 2006; Lai et al. 2007; Finkelstein et al. 
2007; Nilsson et al. 2007, obtained from SED fitting). 

Besides altering the shape of the emitted spectrum, the 
IGM also has an effect on another quantity of much inter- 
est to observers, namely the observed fraction / Q b s of the 
intrinsically emitted number of Lya photons. Since the 
bulk of the emitted Lya photons is due to young stars, 
the total Lya luminosity of a galaxy may be used as a 
proxy for its star formation rate (SFR), one of the main 
quantities characterizing galaxies. Assuming case B re- 
combination, L-^ya can be converted to a total Ha lumi- 
nosity Lhq (through Lh q = L^ ya /8.7; Osterbrock 1989), 
which in turn can be converted to an SFR (through SFR 
= 7.9 x 10" 42 iHa/(erg s" 1 ) Af yr" 1 ; Kennicutt 1998). 

The above conversion factors assume that none of the 
emitted light is lost. If dust is present in the galaxy a 
fraction of the emitted photons will be absorbed (possi- 
bly making the LAE observable in FIR; see Dayal et al. 
2010a). This can be corrected for if the color excess 
E(B — V) is measured, assuming some standard extinc- 
tion curve. However, this assumes that both the Ha and 
Lya radiation are simply reduced by some factor corre- 
sponding to having traveled the same distance through 
the dusty medium. But since the path of Lya photons 
is increased by resonant scattering, this may be far from 
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Fig. 2. — Normalized transmission -F(A) at wavelengths around Lyo, for different redshifts given by the color. In the left panel, the 
vertical axis is linear, while in the right it is logarithmic, emphasizing the transmission at high redshifts. 
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causes a severe suppression of the Lyo line, by z = 3.5 the state of the IGM is not very different. 
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Fig. 5. — Transmission F for z = 3.5 (left), 5.8 (middle), and 6.5 (right), for three different size categories of galaxies; small (dotted), 
intermediate (dashed), and large (solid lines). Although slightly more absorption is seen in the vicinity of smaller galaxies, the transmission 
functions are quite similar for the three size ranges. 
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Fig. 6. — Illustration of the effect of the IGM on the observed Ly« profile emerging from a galaxy at z ~ 3.5. Without taking into account 
the IGM the two peaks are roughly equally high (left panel). However, when the spectrum is transmitted through the IGM characterized 
by the transmission function F(X) (middle panel), the blue peak is dimished, resulting in an observed spectrum with a higher red peak 
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Fig. 7. — Emitted spectra (dotted lines) from nine different simulated galaxies at ~3.5 — ordered after increasing size — and the 
corresponding spectra after being transmitted through the IGM (solid lines with gray regions denoting the 68% confidence intervals). The 
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the truth. Comparing the Lya-inferred SFR with that 
of Ha (or UV continuum), the effect of scattering can 
be constrained, as the quantity SFR(Lya)/SFR(Ha;UV) 
will be an estimate of ,/obs- Using this technique, values 
from a few percent (mostly in the nearby Universe; e.g. 
Hayes et al. 2007; Atek et al. 2008; Hayes et al. 2010) to 
~l/3 at high redshifts (Gronwall et al. 2007) have been 
found. 

Once the radiation enters the IGM, it is also not af- 
fected in the same way, since the IGM is transparent to 
Ha, but not to Lya. Dijkstra et al. (2007) estimated an- 
alytically the fraction of Lya photons that are scattered 
out of the line of sight by the IGM (at z ~ 4.5-6.5) and 
found a mean transmission of /igm ~ 0.1-0.3, depend- 
ing on various assumptions (note that in Dijkstra et al. 
(2007) this fraction is called "T a "). In that study, the 
intrinsic Lya line profile is modeled as a Gaussian, the 



width of which is given by the circular velocity of the 
galaxies, in turn given by their mass, and assuming that 
no dust in the galaxies alters the shape of the line before 
the light enters the IGM. 

Since the transmission function is a non-trivial func- 
tion of wavelength, the exact shape and width of the 
Lya line profile is important. With the Lya RT code 
MoCaLaTA more realistic spectra can be modeled, and 
applying the transmission function found in Section 4.1, 
/igm can be calculated for the sample of simulated galax- 
ies. The transmitted fraction is not a strong function of 
galaxy size; on average, a fraction /igm = 0.77lg;34 of 
the photons escaping the galaxies is transmitted through 
the IGM, at the investigated redshift of z ~ 3.5. The 
results for the individual galaxies is given in Table 2. 

At higher redshifts the metallicity is generally lower, 
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TABLE 2 

Lyo transmission fractions for nine simulated 
galaxies atz~ 3.5 



TABLE 3 

LYa TRANSMISSION FRACTIONS FOR SIX SIMULATED 



Galaxy V c 



/iGM 



S108sc 

S108 

S115sc 

S87 

S115 

K33 

S29 

K15 

S33sc 

Average 



17 
33 
12 
69 
73 
126 
137 
164 
228 



0.97 ±0.02 
0.12 ± 0.02 
0.95 ± 0.03 
0.22 ± 0.01 
0.30 ± 0.05 
0.08 ± 0.02 
0.12 ± 0.03 
0.17 ±0.02 
0.08 ± 0.02 



0.66: 
0.78: 
0.76: 
0.81 
0.80 
0.78 
0.75 
0.79 
0.80; 

0.77^ 



-0.24 
-0.36 



0.35 
0.18 



+0.15 
0.32 



+0.16 
0.34 



+0.17 
0.34 



-0.19 
-0.35 



+0.17 
0.37 



-0.16 
-0.33 



0.64 

0.09 



+0.23 



+0.02 



0.72 



0.18 



0.24 



+0.17 
0.32 

+0.03 
0.07 

+0.06 



0.06 



0.09 



+0.02 
0.03 

+0.03 
0.05 

+0.03 



0.13 



06+ - 02 
u ' UD -0.03 



Note. — Columns are, from left to right: galaxy 
name, circular velocity V c in km s , fraction / csc of 
emitted photons escaping the galaxy (i.e. not absorbed 
by dust), fraction /igm of these transmitted through the 
IGM, and resulting observed fraction / b s = /esc/iGM- 
Uncertainties in / csc represent varying escape fractions 
in different directions, while uncertainties in /igm rep- 
resent variance in the IGM. Uncertainties in / |-> a are cal- 
culated as (o- obs // obB ) 2 = (o- C ac//csc) 2 ±(o-IGM//lGM) 2 - 

leading to less dust and hence larger values of / osc . How- 
ever, the increased neutral fraction of the IGM scatters 
a correspondingly higher number of photons out of the 
line of sight, resulting in a smaller total observed frac- 
tion. Figure 8 and Figure 9 shows the impact of the 
IGM on Lya profiles at z ~ 5.8 and z ~ 6.5, respec- 
tively, and Table 3 and Table 4 summarizes the ob- 
tained fractions. At these rcdshifts, for the six galax- 
ies for which the calculations have been carried out 
on average a fraction /igm(^ = 5.8) = 0.26i"o'ls and 

/igm(z = 6.5) = 0.20lo is °f t ne photons is transmitted 
through the IGM, consistent with what was obtained by 
Dijkstra et al. (2007). 

4.3. Probing the Epoch of Reionization 

We now focus on a different topic, namely the RT in 
the IGM far from the emitting galaxies. Measuring the 
average transmission in a wavelength interval blueward 
of the Lya line for a sample of quasars or other bright 
sources, one can ascertain the average transmission prop- 
erties and hence the physical state of the IGM. The in- 
terval in which the transmission is calculated should be 
large enough to achieve good statistics, but short enough 
that the bluest wavelength does not correspond to a red- 
shift epoch appreciably different from the reddest. Fur- 
thermore, in order to probe the real IGM and not the 
quasar's neighborhood, the upper limit of the wavelength 
range should be taken to have a value somewhat below 
Ac- 

Songaila (2004) measured the IGM transmission in the 
LAFs of a large sample of quasars with redshifts between 
2 and 6.5, in the wavelength interval 1080-1185 A. In Fig- 







GALAXIES AT 


z ~ 5.8 




Galaxy 


V c 


/esc 


/igm 


Jobs 


K15g 


63 


0.93 ±0.02 
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o.33i°ri 5 3 


K15c 
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0.52 ±0.08 
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0.17 ±0.05 


0.16+°- 




Average 






0.26+°;lt 




Note. - 


Same as Table 2, but for z = 5. 
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TABLE 4 

Ly« transmission fractions for six simulated 
galaxies atz~ 6.5 


Galaxy 


V c 


/esc 


/igm 


Jobs 


K15f 


50 


0.96 ±0.03 


24+ ' 18 
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14+ 006 
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0.20i°;°? 
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0.18ig;S« 
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u - uo -0.04 
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Note. — Same as Table 2. but for 



6.5. 



urc 10 we compare the simulated transmitted fractions 
with her sample. As is evident from the figure, an "early" 
reionization, i.e. with z re = 10 (Model 1 and 2), yields a 
slightly too transparent Universe, while a "late" reioniza- 
tion (Model 3) yields a too opaque Universe. Note that 
the log scale makes the z re = 6 data seem much farther 
off than the z rc = 10 data. 

4.4. Convergence test 

In Section 3.4 we stated that the sightlines should ini- 
tiate at a distance ?*o from the centers of the galaxies 
given by their virial radius. Figure 11 displays the cu- 
mulative probability distributions of the distance of the 
last scattering from the galaxy center (calculated with 
MoCaLaTA) for different redshifts and mass ranges, 
demonstrating that in most cases a photon will have ex- 
perienced its last scattering at the order of 1 r v ; r from 
the center of its host galaxy. For increasing redshift, the 
photons tend to escape the galaxies at larger distances 
due to the higher fraction of neutral hydrogen, but the 
change with redshift seems quite slow. Furthermore, at 
a given redshift the distinction between different galactic 
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Fig. 8. — Same as Figure 7, but for z = 5.8. 
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Fig. 9. — Same as Figure 7, but for z = 6.5. 
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Fig. 10. — Comparison of observations (black data points) and 
simulations (colored data points) of the transmitted flux blucward 
of the Lya line as a function of redshift. The three models are 
represented by the colors blue, green, and red for model 1, 2, and 
3, respectively. To highlight the significance of the improved UV 
RT, we show both the transmission in the "original" simulation 
with the "pseudo"-RT (crosses) and with the improved UV RT 
(squares). For details on the observations sec Songaila (2004), 
from where the data are kindly supplied. 

size ranges appears insignificant (as long as ro is mea- 
sured in terms of r v j r ). 

Also shown in Figure 11 are the transmission curves 
for sightlines initiating at various distances ro from the 
centers of the galaxies. For very small values of ro, a 
significantly lower fraction is transmitted due to the high 
density of neutral gas. Around r ~ r v i r the change in 
F(X) becomes slow, converging to F(X) exhibiting no 
dip for ro — » oo. 

In summary, scaling ro to the virial radius r v ; r of the 
galaxies allows us to use the same value of ro = 1.5r v ; r 
for all sightlines. 

5. DISCUSSION 

5.1. The origin of the dip 

In general, the effect of the IGM — even at relatively 
low redshifts — is to reduce the blue peak of the Lya 
line profile. At z ~ 3.5, the effect is not strong enough to 
fully explain the oft-observed asymmetry, but at z > 5, 
almost all radiation blucward of the line center is lost in 
the IGM. 

Figures 12 and 13 display for three different epochs the 
velocity field and the density of neutral hydrogen, respec- 
tively, of the IGM associated with the galaxies, taken as 
an average of all galaxies in the samples, and in all di- 
rections. At all redshifts, on average the IGM close 
to the galaxies recedes at a somewhat slower rate than 
that given by the pure Hubble flow, slowly converging to- 
ward the average expansion rate of the Universe. At the 
higher redshifts, the IGM in the immediate proximity of 
the galaxies is characterized by higher velocities, due to 
outflows generated by starbursts. At z = 3.5, however, 
this effect is overcome by the accretion of gas. 

Inspecting the dips in Figure 2, the minima are seen to 
be located at roughly 50 km s _1 , almost independently 



of the redshift but becoming slightly broader with in- 
creasing z. At z = 3.5, the dip extends all the way out 
to 300-400 km s _1 . This corresponds to the central ab- 
sorption being caused by the IGM within ~150 kpc, and 
the wings of the absorption by the IGM within ~1 Mpc. 
As shown in Figure 13, at z = 3.5 the density of neutral 
hydrogen decreases monotonically with distance from the 
source, and within ^150 kpc Uhi is substantially higher 
than the cosmic mean. Further away, the density is close 
to the mean density of the Universe. However, as seen in 
Figure 12 the recession velocity of the gas continues to 
lie below that of the average, Universal expansion rate, 
and in fact does so until approximately 1 Mpc from the 
source. 10 

Thus, the cause of the suppression of the blue wing of 
the Lya line may, at wavelengths close to the line center 
(AA ~ 1/2 A), be attributed to an increased density 
of neutral hydrogen close to the galaxies, while farther 
away from the line center (AA < 1—1.5 A), to a retarded 
Hubble flow. 

Zheng et al. (2010a) perform Lya RT in the IGM at 
z = 5.7 through full Monte Carlo simulations, and make 
a detailed comparison with the observed line profile ob- 
tained from simply multiplying the intrinsic line by e _T , 
as has been done in previous models (e.g. Iliev et al. 
2008). They conclude that neglecting scattering effects 
severely underestimates the transmitted fraction. While 
it is certainly true that treating scattering processes as 
absorption inside the galaxies is only a crude approxima- 
tion, once the probability of photons being scattered into 
the line of sight becomes sufficiently small, this approach 
is quite valid. In their analysis, Zheng et al. (2010a) start 
their photons in the center of the galaxies, which are re- 
solved only by a few cells (their dx being ^28 kpc in 
physical coordinates and their fiducial galaxy having r v i r 
= 26 kpc). Since the side length of our smallest cells 
are more than 400 times smaller than the resolution of 
Zheng et al. (2010a), we are able to resolve the galax- 
ies and their surroundings in great detail, and we are 
hence able to determine the distance at which the e~ T 
model becomes realistic. Moreover, when coupling the 
IGM RT with Lya profiles, we use the realistically cal- 
culated profile, whereas Zheng et al. (2010a) use a Gaus- 
sian set by the galaxies' halo masses. Their line widths 

Omit = 32A/ 1 1 g 3 km s" 1 , where Mio is the halo mass 
divided by lO lo /i _1 M , thus neglect broadening by scat- 
tering. This makes them much smaller than ours, which 
are typically several hundred km s _1 . Note, however, 
that our relatively small cosmological volume and the 
fact that ionizing UV RT is performed as a post-process 
rather than on the fly may make our density field less 
accurate than that of Zheng et al. (2010a). 

The photons that are scattered out of the line of sight 
are of course not lost, but rather become part of a dif- 
fuse Lya background. Since more scatterings take place 
in the vicinity of galaxies, LAEs tend to be surrounded 
by a low-surface brightness halo, making them look more 
extended on the sky when comparing to continuum bands 

10 Since the cosmological volume is several Mpc across, except 
for the galaxies lying close to the edge, the absorption takes place 
before the sightlines "bounce", so in general there is no risk of 
a sightline going through the same region of space before having 
escaped the zone causing the dip. 
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Fig. 11. — Left: Cumulative probability distribution of the distance r from the center of a galaxy at which the last scattering takes place 
(calculated with the Lya RT code MoCaLaTA (Laursen et al. 2009a)). Thin solid lines represent individual galaxies at redshift 2.5 (blue), 
3.5 (green), and 6.5 (red), while thick solid colored lines are the average of these. Also shown, in black, are the average of three different 
size ranges at z = 3.5; small (dotted), intermediate (dashed), and large galaxies (solid). 

Right: Resulting transmission function -F(A) for sightlines originating at various distances ro/r v i r from galactic centers, increasing in steps 
of 0.1r v ; r and ranging from ro = 0.1r v ; r to ro = 5r v ; r . The three different redshifts are shown in shades of blue, green, and red, for z = 2.5, 
3.5, and 6.5, respectively, and the lines go from light shades for low values of ro to dark shades for high values of ro- The transmission 
functions corresponding to ro = 1.5r v ; r are shown in black dotted lines. All results arc for Model 1. 
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Fig. 12. — Average recession velocity i>b u ik of the IGM as a function of proper distance d from the centers of the galaxies in Model 1 
(solid black lines, with gray regions indicating the 68% confidence intervals). At all redshift, the expansion is retarded compared to the 
pure Hubble flow (dashed line) out to a distance of several comoving Mpc. At high redshifts, however, very close to the galaxies outflows 
generate higher recession velocities. 
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Fig. 13. — Average density rijjj of neutral hydrogen as a function of proper distance d from the centers of the galaxies in Model 1 
(solid black lines, with gray regions indicating the 68% confidence intervals). While in general the density decreases with distance, at high 
redshifts ionizing radiation reduces njj : in the immediate surroundings of at least some of the galaxies, as seen by the small dip at ~50 
kpc in the lower part of the gray zone. 
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on the r ~ 10-100 kpc scale. This has been confirmed 
observationally (e.g. M0ller & Warren 1998; Fynbo et al. 
2001, 2003, but see Bond et al. (2010)), and studied the- 
oretically/numerically (e.g. Laursen & Sommer-Larsen 
2007; Barnes & Haehnelt 2010; Zheng et al. 2010c). As 
argued above, part of the cause of the line suppression 
(mostly in the wings) is due to the IGM up to ~1 Mpc 
of the source. At this distance, the surface brightness 
is much lower than observational thresholds, but could 
be detected by stacking images of LAEs (Zheng et al. 
2010c) 

As seen in e.g. Figure 3, a large scatter between in- 
dividual sightlincs exists, reflecting the generally quite 
inhomogeneous IGM. Consequently, one cannot simply 
use the calculated transmission function for deconvolv- 
ing observed Lya lines to obtain the intrinsic line profiles, 
other than in a statistical sense. With a large sample, 
however, more accurate statistics on Lya profiles could 
be obtained. Calculating the transmission functions as 
an average of all directions, as we have done in this work, 
assumes that observed galaxies are randomly oriented in 
space, i.e. that there is no selection effects making more 
or less luminous directions pointing toward the observer. 
For LAEs clustered in, e.g., filaments, the effect of the 
retarded Hubble flow may be enhanced perpendicular to 
the filament, making the galaxies more luminous if ob- 
served along a filament than perpendicular to it (see Fig- 
ure 9 in Zheng et al. 2010b). 

5.2. Galactic outflows 

As discussed in the introduction, at high redshifts 
many galaxies are still in the process of forming and are 
expected to be accreting gas. In principle, this should re- 
sult in a blueshifted Lya profile, but this is rarely seen. 
Evidently, IGM absorption is unable to always be the 
cause of this missing blue peak. On larger scales, mass 
is observed to be conveyed through large streams of gas; 
the cosmic filaments. Although this has been seen only 
tentatively on galactic scales, galaxy formation may be 
expected to occur in a similar fashion. Indeed, numerical 
simulations confirm this scenario (e.g. Dekel et al. 2009; 
Goerdt et al. 2010). Very recently, Cresci et al. (2010) 
reported on an "inverted" metallicity gradient in three 
z ~ 3 galaxies, which they interpret as being due to the 
central gas having been diluted by the accretion of pri- 
mordial gas. In contrast to gas accreting through a few 
narrow streams, outflows are more or less isotropic. Even 
bipolar outflows have a rather large opening angle, and 
thus the probability of a sightlinc towards a galaxy pass- 
ing through outflowing gas is larger than passing through 
infalling gas. Thus, it may be expected that more obser- 
vations are obtained of Lya profiles lacking the blue peak 
than lacking the red peak. 

The fact that starbursts are needed to generate large 
outflows also imposes a bias on the observations; even 
if outflows happen only during relatively short phases in 
the early life of a galaxy (for LBGs, Ferrara & Ricotti 
(2006) found that a typical starburst phase lasts only 
about 30±5 Myr), such galaxies are more likely to appear 
in surveys, simply because they are more luminous than 
those with small SFRs. 

5.3. Transmitted fraction of Lya photons 



Even though absorption in the IGM does not alter the 
line shape drastically at an intermediate redshift of 3.5, 
it reduces the intensity by roughly one-fourth, as seen 
from Table 2. This fraction is not a strong function of 
the size of a galaxy, but since the spectra emerging from 
larger galaxies tend to be broader than those of smaller 
galaxies, a comparatively larger part of the small galaxy 
spectra will fall in the wavelength region characterized by 
the dip in the transmission function. Hence, on average 
the IGM will transmit a larger fraction of the radiation 
escaping larger galaxies. 

One- fourth is not a lot, but since it is preferentially 
blue photons which are lost, the spectrum may become 
quite skewed when traveling through the IGM. The lost 
fraction /igm is in addition to what is lost internally in 
the galaxies due to the presence of dust. As mentioned 
in Section 4.2, dust tends to make the line profile more 
narrow. For galaxies with no dust, the lines can be very 
broad. In this case, /igm will be slightly higher, since for 
broad lines a relatively smaller part of the spectrum falls 
on the dip seen in F( A) . Note however that in general dis- 
tinguishing between dust absorption and IGM effects will 
be rather tricky, if not impossible (see also Dayal et al. 
2010a). 

As expected, at higher redshifts the IGM is more 
opaque to Lya photon; at z — 5.8 and z = 6.5 the inten- 
sity was found to be reduced by approximately 75% and 
80%, respectively. At z = 6.5 the blue wing is completely 
lost. In most cases this is true also at z = 5.8. However, 
as is seen from the gray area representing the 68% con- 
fidence interval, in some cases an appreciable fraction of 
the blue wing can make it through the IGM. 

5.4. The significance of dust 

The IGM calculations performed in this work are all 
neglecting the effect of dust. However, dust is mostly 
(although not completely) confined to the galaxies, es- 
pecially the central parts, and since by far the greatest 
part of the distance covered by a given sightline is in the 
hot and tenuous IGM, one may expect this to be a fair 
approximation. 

This anticipation was confirmed following the dust 
model of Laursen et al. (2009b). In short, dust density 
scales with gas density and metallicity, but is reduced in 
regions where hydrogen is ionized to simulate the effect of 
dust destruction processes. The factor tihic(A) in Equa- 
tion (2) is then replaced by nHic(A) + rid<Td(A), where 
rid and era (A) are the density and cross section of dust, 
respectively. The resulting decrease in transmission is at 
the 10~ 4 to 10~ 3 level. 

5.5. "Early" vs. "late" reionization 

As is evident from Figure 10, although none of the 
models really fit, the early and late reionization bracket 
the observations. Model 1, with erg = 0.74 provides a 
somewhat better fit to the observations than Model 2 
with og = 0-9- Due to the more clumpy structure of the 
latter, galaxies tend to form earlier, rendering the IGM 
more free of gas and thus resulting in a slightly more 
transparent Universe. However, this should not be taken 
to mean that the lower value of a% is more realistic, since 
a higher £7§ could be accounted for by a (slightly) smaller 
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Investigating the LAF is not the only way of prob- 
ing the EoR. The WMAP satellite detects the effects 
of the Thomson scattering of cosmic microwave back- 
ground (CMB) radiation on free electrons. Since the 
number density n e of free electrons increases as the Uni- 
verse gets reionized, a signature of the EoR can be ob- 
tained by measuring the total optical depth r e to Thom- 
son scattering. The latest results (Jarosik et al. 2010) 
give r e = 0.088 ± 0.015, implying for their model of the 
EoR that the Universe was reionized at z ro i on ~ 10-11. 

The optical depth to Thomson scattering in a short 
interval dt of time is 

dr = n e (z) ot cdt 

= n e (z)a T c 1 dz, (3) 
(l + z)H(z) 

where ctt is the Thomson scattering cross section. That 
is, given an electron density history, the resulting total 
r e can be calculated by integrating Equation (3). 

Figure 14 displays the average comoving number den- 
sity n e of electrons as a function of redshift, as well as the 
corresponding total optical depth r e of electrons. The re- 
sults for Model 2 and Model 3 are shown; those for Model 
1 lie very close to those of Model 2. Going from larger 
toward smaller rcdshifts, both models are seen to be char- 
acterized by a roughly constant and very low density of 
electrons before reionization, then a rapid increase not 
at, but shortly after z re , and finally a slow decrease, due 
to subsequent gas cool-out and resulting galaxy forma- 
tion. The rapid increase marking the EoR lasts approxi- 
mately 100 Myr. The difference in ionization history for 
the pseudo- and the realistic UV RT schemes is not criti- 
cal, although in the latter the creation of ionized bubbles 
around stellar sources causes a slightly earlier EoR than 
the pseudo- RT is able to, especially in Model 3 (z IC = 6). 

Also shown in Figure 14 is the corresponding r e history 
that would prevail in the hypothetical case of an instant 
reionization, where the Universe is fully neutral before, 
and fully ionized after, some redshift z re i on fixed to make 
the total T e match the value measured by WMAP. In this 
model, no gas is assumed to be locked up in stars, and 
n e is thus given by 

f for Z > Z rG i on 

ne{z) = i fOT * < 

where ip = 0.76 is the mass fraction of hydrogen, 
Qb = 0.046 is the baryonic energy density parameter 
(Jarosik et al. 2010), p c is the critical density of the Uni- 
verse, and mn is the mass of the hydrogen atom. 

As seen from the figure, not even the early reioniza- 
tion model is able to reproduce the optical depth probed 
by WMAP. However, for a model with an even earlier 
EoR, the transmission T of the IGM would be even far- 
ther off the observational data, as seen from Figure 10. 
In general, there seems to be a significant disagreement 
between the WMAP results and the QSO results con- 
cerning the redshift of the EoR. Many authors have tried 
to resolve this apparent discrepancy, e.g. Wyithe & Loeb 
(2003) and Cen (2003) who considered a double reioniza- 
tion, first at z ~ 15-16 and later at z ~ 6. 

Notice nevertheless that the total r e of the z re = 10 
model is not entirely inconsistent with the WMAP- 
inferred value, being roughly 2a away. Furthermore, 



since only the total optical depth is measured by WMAP, 
the exact history of the EoR is obviously less certain. 
Generally, either an instant reionization must be as- 
sumed, or perhaps a two-step function to make the EoR 
slightly extended, and possibly with an additionally step 
at z ~ 3 to account for helium reionization. 

The above results are at odds with the interpretation 
of the CMB EE polarization maps as showing that the 
reionization of the Universe was complete at z = 10.5 
(Jarosik et al. 2010), regardless of which model is as- 
sumed for the EoR. However, many factors enter the con- 
version of the polarization maps into an optical depth, 
and the model of this EoR may be too simplistic. Com- 
pared to WMAP, the recently launched Planck 11 satel- 
lite has a much higher resolution and sensitivity; when 
Planck data become available, these issues may be solved 
as much tighter constraints can be put on parameters like 
t € (Galli et al. 2010). 

6. SUMMARY 

Simulating sightlines through cosmological hydrosim- 
ulations, we have investigated the transmission through 
the IGM of light in the vicinity of the Lya line. The high 
resolution of the cosmological simulations combined with 
the adaptive gridding for the RT allows us to probe the 
vclcity field around the galaxies in great detail. While 
earlier studies (Zheng et al. 2010a) have shown that this 
approach of simply multiplying an e~ T factor on the 
intrinsically emitted Lya line is a poor approximation 
for the observed line profile and transmitted fraction, 
we have argued that when the circumgalactic environs 
are sufficiently resolved and combined with realistically 
calculated intrinsic Lya lines, this approach should be 
valid. Special emphasis was put on how Lya line profiles 
emerging from high-redshift galaxies are reshaped by the 
surrounding IGM. In general a larger fraction of the blue 
side of the line center Ao is lost as one moves toward 
higher rcdshifts. At z > 5, almost all of the light blue- 
ward of Ao is lost, scattered out of the line of sight by the 
high neutral fraction of hydrogen. However, even at rela- 
tively low redshift more absorption takes place just blue- 
ward Ao- A transmission function F(X) was calculated, 
giving the fraction of light that is transmitted trough the 
IGM at various epochs. At all redshifts where some of 
the light blueward of Ao is transmitted (i.e. at redshifts 
below ~5.5), a significant dip with a width of the order 
1 A is seen. The origin of this dip is a combination of 
an increased density of neutral hydrogen and a retarded 
Hubble flow in the vicinity of the galaxies. 

This extra absorption may in some cases be the reason 
that the blue peak of an otherwise double-peaked Lya 
profile is severely reduced, or lacking. Nevertheless, it is 
not sufficient to be the full explanation of the numerous 
observations of z ~ 2-4 Lya profiles showing only the 
red peak. The outflow scenario still seems a credible 
interpretation. 

Combining the inferred transmission functions with 
simulated line profiles, the fraction of Lya photons that 
arc transmitted through the IGM at z ~ 3.5, 5.8, and 
6.5 was computed, and was found to be /igm(z = 3.5) = 
0-77i°'iI, fi GM (z = 5.8) = 0.26±g;il, and f IGM (z = 

11 http:/ / www.rssd.esa.int /Planck/ 
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Fig. 14. — Left: Volume averaged, comoving electron density n e /(l + .z) 3 as a function of redshift z for Model 2 (green) and 3 (red) which 
have z TE = 10 and z re = 6, respectively. Both models arc characterized by a quite sharp increase in n e shortly after the onset of the UVB. 
For Model 2 the EoR is seen to take place around a ~ 8.5, while for Model the EoR lies at z ~ 5.5. 

Right: Integrated optical depth r e of electrons as a function of redshift z for the two models. The dark gray line with the associated lighter 
gray 68% and 95% confidence intervals indicates the T e history required to reach the optical depth measured by WMAP, if an instant 
rcionization is assumed. As expected, none of the models are able to reach the 0.088 ± 0.015 inferred from the WMAP results, although 
the z Ie = 10 model is "only" ~2cr away. Since the n e data do not extend all the way to z = 0, a fiducial value of n e (z = 0) = 1.5 X 10 — 7 
cm -1 has been used. The exact value is not very imporant, since the proper density at low redshift is very small. At low redshift, the 
model curves lie slightly below the theoretical curve. This is due to a combination of the models including helium ionization, releasing 
more electrons, and star formation and gas cool-out, removing free electrons. 



6.5) = 0.20to'ig, respectively. This is in addition to what 
is lost internally in the galaxies due to dust. The stan- 
dard deviations were found to be dominated by sightlinc- 
to-sightline variations rather than galaxy-to-galaxy vari- 
ations. 

Considering the average fraction of light far from the 
line on the blue side transmitted through the IGM, and 
comparing to the comprehensive set of observations of 
the LAF by Songaila (2004), we constrain the EoR to 
have initiated between z rc = 10 and z Te = 6, correspond- 
ing the having ionized a significant fraction of the Uni- 
verse around z ~ 5.5 and z ~ 8.5, respectively. Even 
though the "early" models of z rc = 10 produces a slightly 
too transparent Universe when comparing to the LAF, 
the optical depth of electrons is too low when comparing 
to the observations of the CMB by the WMAP satellite, 



possibly indicating a too simplistic interpretation of the 
CMB polarization. 
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